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We consider the nonlinear Dirac equations (NLDE's) in 1+1 dimension with scalar-scalar self 

2 - i 

interaction ^^py^^) in the presence of various external electromagnetic fields. We find exact 
solutions for special external fields and we study the behavior of solitary wave solutions to the NLDE 
in the presence of a wide variety of fields in a variational approximation depending on collective 
coordinates which allows the position, width and phase of these waves to vary in time. We find that 
in this approximation the position q(t) of the center of the solitary wave obeys the usual behavior 
of a relativistic point particle in an external field. For time independent external fields we find that 
the energy of the solitary wave is conserved but not the momentum which becomes a function of 
time. We postulate that similar to the nonlinear Schrodinger equation (NLSE) that a sufficient 
dynamical condition for instability to arise is that dP(t)/dq(t) < 0. Here P(t) is the momentum 
of the solitary wave, and q is the velocity of the center of the wave. We investigate the accuracy 
of our variational approximation using numerical simulations of the NLDE and find that when the 
forcing term is small and we are in a regime where the solitary wave is stable, that the behavior of 
the solutions of the collective coordinate equations agrees very well with the numerical simulations. 

PACS numbers: 05.45.Yv, 03.70. +k, 11.25.Kc 



I. INTRODUCTION 



Classical solutions of nonlinear field equations have a long history as a model of extended particles [2 [2]. In 1970, 
Soler [2] proposed that the self-interacting 4-Fermi theory was an interesting model for extended fermions. Later, 
Strauss and Vasquez [3J were able to study the stability of this model under dilatation and found the domain of 
stability for the Soler solutions. Solitary waves in the 1+1 dimensional nonlinear Dirac equation (NLDE) have been 
studied [4] [5] in the past in the case of massive Gross-Neveu [6] (with N = 1, i.e. just one localized fermion) and 
massive Thirring [7] models). In those studies it was found that these equations have solitary wave solutions for 
both scalar-scalar (S-S) and vector-vector (V-V) interactions. The interaction between solitary waves of different 
initial charge was studied in detail [5] for the S-S case in the work of Alvarez and Carreras by Lorentz boosting the 
static solutions and allowing them to scatter. Recently we extended the solutions previously found to a more general 
interaction of the form ^rpj( , & v I') K+1 [9]. For the nonrelativistic limit of the NLDE, namely the nonlinear Schrodinger 
equation (NLSE), there have been recent studies of the behavior of the NLSE in external fields. Using a collective 
coordinate theory, the authors found [TU] [TT], [T3] [13] that a sufficient dynamical condition for instability to arise 
is that dp(t)/dv < 0. Here p(t) is the normalized canonical momentum p{t) = = J dx^*(x, t)^(x, t) is 

the mass and q(t) = v(t) is the velocity of the solitary wave. 

One of the points we will investigate in the paper is whether this dynamical stability criterion is also valid for the 
NLDE. There has been recent interest in the stability of NLDE with higher-order nonlincarity [H] • Comech (private 
communication) has been able to prove that for k = 1, the Vakhitov-Kolokolov [15] criterion guarantees linear stability 
in the non-relativistic regime of the NLDE equation for solutions of the form (in the rest frame) *5>(x,t) = il)(x)e~' iu)t 
where uj is less than but approximately equal to the mass paramater m in the Dirac equation. He was also able to 
show linear instabilty in the same nonrelativistic regime for k > 3. This is the first rigorous result for the Dirac 
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equation, but it only applies in the nonrelativistic regime. Here we want to understand if we can determine in the 
relativistic regime for what values of w do the solitary waves become unstable, with and without forcing terms even 
when they arc stable in the non-relativistic regime. 

This paper is organized as follows: In Sec. II, we review the known exact solutions and discuss the conservation 
laws that govern their behavior. In Sec. Ill we extend Bogolubsky's discussion [TB] of the stability of solitary waves 
to arbitrary nonlinearity parameter n. In Sec. IV we consider the NLDE in external electromagnetic fields, find 
particular exact solitary wave solutions, and discuss the stability of these solutions. In Sec. V we introduce our 
variational method based on using the exact solution wave functions for the unforced problem, while having the 
position, width parameter and phase become collective coordinates depending on time. We write the relativistic 
equations for these collective coordinates which are similar to point particle relativistic dynamical equations. In 
Sec. VI we postulate our stability criteria for an arbitrary external potential. In Sec. VII we examine and solve 
the collective coordinate (CC) equations for three types of potentials-a ramp potential, a harmonic potential and a 
spatially periodic potential. We also compare the solution to the CC equations to the numerical simulation of the 
NLDE equation. In the appendix we discuss identities that are obeyed by the solutions in the rest frame. 

II. REVIEW OF EXACT SOLUTIONS TO THE NLDE 

In this section we review the exact solutions to the NLDE, using the notation of [5]. We are interested in solitary 
wave solution of the NLDE given by 

-m)* + g 2 (**) K * = 0. (2.1) 
These equations can be derived in a standard fashion from the Lagrangian 

L = (-} [#7^* - dj>^] ~ to** + (**) re+1 . (2.2) 

\2 J K + 1 

For solitary wave solutions, the field * goes to zero at infinity. It is sufficient to go into the rest frame, since the 
theory is Lorentz invariant and the moving solution can be obtained by a Lorentz boost. In the rest frame we consider 
solutions of the form 

¥(x,t) = e~ iut i){x). (2.3) 

We are interested in bound state solutions that correspond to positive energy oj > and which have energies in the 
rest frame less than the mass parameter m, i.e. u> < to. In our previous paper [9], we chose the representation 70 = (73, 
«7i = G\. Here instead, to make contact with the numerical simulations paper of Alvarez and Carreras [8] we choose 
instead 



Defining A, B via: 



7° = ^3; 7 1 = ^2- (2.4) 



we obtain the following equations for A and B. 

HA 

+ (to + u)B - g 2 (A 2 - B 2 ) K B = , 

ax 

^ + (m - u)A - g 2 (A 2 - B 2 ) K A = . 
ax 

(2.6) 

A first integral of these equations can be obtained by realizing that from energy-momentum conservation we have 
Thus for stationary solutions 

T 10 = const. T 11 = const. (2.8) 



Now using (2.3 1 we obtain 



T 11 = ui/rij) — miptp + Lj, 



where 



9 



K + l 



WO 



K+l 



For solitary wave solutions vanishing at infinity the constant is zero and we get the useful first integral: 

T 11 = Lo^ip - mtfip + L T = 0. 



Multiplying the equation of motion on the left by ^ and using (2.3) we have that: 



(k + l)Lj = —cjip't/j + mipip — ipi^f 1 d\ip . 



Therefore we can rewrite T 11 — as 



ujnip'ij) — rriKipip ~ V'*7 1 ^iV' = 0. 



For the Hamiltonian density we obtain 



n = T m = - [*7 1 9 a; * - d^j 1 ^] + m** -L I = h 1 + h 2 - h 3 . 



Each of the hi are positive definite. From Eq. (2.111 and( 2.12) one has the relationship: 



From this we have 



and in particular for k = 1, H = rmjjip. 
In terms of R, one has 



kLi — —ifii , yd x ip. 



h 3 = -h± 

K 



d8 

dx 



This leads to the simple differential equation for 9 for solitary waves 

d0 



dx 



The solution is 



where 



Thus we have 



-oj k + m K cos 20; where lo k = k cj; m K — Km. 



0(x) — tan 1 (atanh/3 K a;), 





\ 1/2 


/ m — 






+ UJ K J 


\m + 





1/2 



, p K = {mi - uiyi*. 



sin t/(x) 



a tanh j3 K x (m — lj) sinh /3 K x 



cos 2 0{x) = 



1 + a 2 tanh (3 K x m cosh 2f3 K x + lj 
(m + lj) cosh 2 ft K x 



1 



1 + a 2 tanh 2 (3 K x m cosh 2f3 K x + lj 



(2.9) 

(2.10) 

(2.11) 
(2.12) 
(2.13) 

(2.14) 
(2.15) 
(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



where we have used the identities: 



1 + a 2 tanh 2 fi^x — 



1 — a 2 tanh fi^x — 



m cosh 2/3^2; + u) \ 

sech pkX 

m + u) 

w cosh 2/3^3; + to 

m + to 



sech /3kX 



We have from Eq. (2.101 and Eq. ( 2.11) 



loR 2 - mR 2 cos 29 + — — (i? 2 cos 26 



. K+l 



Thus 



Now we have 



(k + 1) (m cos 29 — uj) 
g 2 (cos26) K + 1 



1/k 



de = pi 

dx uj k + m K cosh 2(3 K x 



—uj k + m K cos 29, 



where /3 K = yj to 2 — uj 2 = ny/m 2 — uj 2 , so that 

to k + uj k cosh 2(3 K x m + uj cosh 2(3 K x 
uj k + to k cosh 2(3 K x uj + m cosh 2(3 K x 



cos 26» 



One important expression is 



Using this we get 



\ cos 29 — uj = 



o / w + to cosh 2/3 K a: 
ri = 



to + uj cosh 2/3^2; 



k 2 (cj + to cosh 2(3 K x) 



g 2 K 2 (m + uj cosh 2j3 K x) 



i/k 



Using the identities in Eq. (2.22), we obtain the alternate expression 



R z 



1 + a 2 tanh /3 K x 
1 — a 2 tanh /3 K a: 



sech 2 /3 K a;(K + l)/3 2 



-i 1/k 



' 2 K 2 (m + w)(l — a 2 tanh /3 K ar) 



In particular for n = I 



and 



9 2(m-w) (1 + a 2 tanh /3a;) , ,„ 

i? 2 = — = — ^— ^sech 2 ^a; 

5 2 (1 - a 2 tanh 2 /3a;) 2 

2 (to 2 — w 2 ) (to + w) cosh 2 /3a: 



A 2 = R 2 cos 2 6» = 
B 2 = R 2 sin 2 9 = 



g 2 (to + uj cosh 2/3x) 2 
2 (to 2 — cj 2 )(to — w) sinh 2 /3a: 
g 2 (to + u; cosh 2/3 x) 2 



We can also write 



cos( 



(to + uj) cosh 2 (K/3a;) 
w + to cosh(2«/?x) 



sin 6* 



(m — uj) sinh 2 (K/3a:) 
+ TOCOsh(2At/3a;) 
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so that 



A = 



B 



' (m + u) cosh 2 (k/3x) 
m + lu cosh(2n (3 x) 

' (m — w) sinh 2 («;/3a;) 
m + w cos1i(2k/3x) 



( K +l)/3 2 



5 2 (m + w cosh(2«;/3a;)) 
(k + 1)/3 2 



g 2 (m + wcosh(2/c/3x)) 



Because of Lorentz invariance we can find the solution in a frame moving with velocity v with respect to 
frame. The Lorentz boost is given in terms of the rapidity variable r\ as follows (here c = 1): 



v = tanhfy; 7 = 



— cosh -q: sinhr/ = 



In the moving frame, the transformation law for spinors tells us that: 



, . _ / cosh(ry/2) sinh(??/2) \ / *°[ 7 ( x — vt),j(t — vx) 
^ ' ; ~ I sinh(7j/2) cosh(r//2 J\ V%[j{x - vt), -y(t - vx) 



since 



cos 



h( V /2) = v/(l+7)/2; sinhfa/2) = y/fr - l)/2, 



This in component form: 



where 



*i(a;,t) = (cosh(77/2)^(x') +isinh(r7/2)B(x'))e-^*', 
$ 2 (x,t) = {smh(r]/2)A(x')+icosh(r]/2)B(x'))e- iwt ' , 

x' = j(x — vt); t' = -f(t — vx). 



(2.33) 
the rest 

(2.34) 

(2.35) 
(2.36) 

(2.37) 
(2.38) 



Note that cosh 2 (77/2) + sinh 2 (?7/2) = cosh 77 = 7. 



A. Conservation Laws of the NLDE 



The Lagrangian is invariant under the transformation of phase ^ — > e lA \I/, which by Noether's theorem leads to 
the conserved current: 



d^(x)=0; j" = *y*. 



This leads to charge conservation: 



which for the solitary wave solution leads to 

Q = j dx(A 2 + B 2 ) 



1 

«V3 



where 



Ua 2 )= f^y— 



(k + 1)P 2 
g 2 (m + lu) 



1 + a 2 y 2 



1/k 



I K (a 2 ). 



= B 
+a 2 B 



1 1 

2' K 
3 1 

2' K 



j/2)(«-l)/«[l _ 0,2^2] 
111 1 9 

2 ^ 1 + ^2'2 + ^ 



.- , . 13 3 1 2 

k 2 2 K 



(2.39) 



(2.40) 



(2.41) 



(2.42) 
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and a 



2 m— u 



2F1 is a hypergeometric function and B(x, k) denotes the beta function. 



We also have energy- momentum conservation Eq. (2.7) leading to conservation of energy: 



and momentum conservation 



E = I T QQ dx, 



T 01 dx. 



(2.43) 



(2.44) 



Because of Lorentz invariance it is sufficient to calculate the energy momentum tensor in the comoving frame v = 0. 
The energy momentum tensor in an arbitrary frame is then given by 



where 



^ _ f cosh 77 sinhry 
a \ sinh?7 cosh 77 



In the rest frame of the solitary wave, for the unperturbed system one has that 

1' 



where 



Too = hi 1 



K /3 2 



dx m + u> cosh(2«;/3a;) \_g 2 (m + to cosh(2«/?x)) 



(k + 1)/3 2 



1/k 



(2.45) 



(2.46) 



(2.47) 



(2.48) 



h% = rntpijj = m(A 2 — B 2 ) = m 



( K + l)/3 2 



_g 2 {m + uj cosh(2/«/3x)) 
Integrating in the rest frame, we get for the rest frame energy 

1 



1/k 



Eo = Hi 1 - 



H2, 



where 



H\ = / dxhi 



(3 



m + uj 



1/k 



g 2 (m + uj) 
113 1 



7 M^1 + H 2F1 (1 + -, J,^ + -;« 2 

2 k / \ k22k 



(2.49) 



(2.50) 



(2.51) 



Ho 



dxho 



K/3 



(«+l)/3 2 
g 2 (m + w) 



2 iV* 



, 1 1\ /111 1 9 

2 k ) \ K 2 2 k 



(2.52) 



Since in the rest frame for stationary solutions T 11 = T 01 = 0, the energy of the solitary wave in the moving frame is 
just 



E = E cosh 77 = jE ; P = E sinh r), 



(2.53) 



so that the norm E 2 — P 2 = E 2 ~ M$. 
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In particular, for k — 1 and m — 1, we have that 



with Q and w related as follows: 



We also have 



M = ^-sinh- 1 ^, (2.54) 
9 Q 2 



Q=^-\ f3=Vl-0J 2 . (2.55) 
9 u 



Hi = 



2(7r 



^-2tanh- 1 ( v /^f 
9 2 



4tanh _1 ( x f^£ 

H 2 = (2.56) 

g 2 

III. STABILITY OF EXACT SOLUTIONS 
A. Stability to changes in the frequency at fixed charge 

Bogolubsky [TH] suggested that the stability could be ascertained by looking at variations of the wave function, 
keeping the charge fixed and seeing if the solution was a minimum (stable to that variation) or maximum (unstable 
to that variation) of the Hamiltonian as a function of the parameter uj. This principle has been very useful in the 
past to determining the stability of scalar wave equations that are Hamiltonian dynamical systems. If the variation 
decreases the energy it turned out that the solitons were unstable. Since in higher dimensions there are many degrees 
of freedom for perturbing the system, this criterion is a sufficient condition for instability. For the Dirac case we have 
found from our numerical simulations that this criterion does not determine the critical uj except when k = 1 [17] . 
the case originally studied by Bogolubsky [16] . Assuming we know the wave function at the value of uj corresponding 
to a fixed charge Q, if we change the parametric dependence on uj this also changes the charge. This can be corrected 
by assuming that the new wave function has a new normalization that corrects for this. That is if we parametrize a 
rest frame solitary wave solution of the NLDE which has a charge Q[uj] by 

ip s (x,t)= Xs {x,u)e- iwt , (3.1) 
then we choose our slightly changed wave function to be 

ip[x,t,w ,uj\ = -^==Xs(x,uj )e 

^fiuj'^Xs^uj'y^'K (3.2) 

Then the wave function ip[x, t, uj' , uj] has the same charge as ij>[x, t, uj]. Inserting this wave function into the Hamiltonian 
we get a new Hamiltonian H p depending on both uj',uj. As a function of uj' this new Hamiltonian is stationary as 
a function of uj' at the value uj' — uj. The criterion Bogolubsky proposed [16] is that the solitary wave is stable 
(unstable) according to whether this new Hamiltonian has a minimum (maximum) at uj' — uj. What we will find for 
K = 1 is that there is a critical value of uj (determined by the coupling g and Q) below which the solitary wave is 
unstable, and this result is borne out by numerical simulations which we will present below. However, we will present 
in another paper numerical simulations at arbitrary k which suggest that this method gives a sufficient condition for 
instability [17 . The probe Hamiltonian has the form: 

H p [uj\uj]=H 1 [uj'] (f(oj',u) 2 - -J{uj',ujf^ +1 ^ +H 2 [uj']f(uj',uj) 2 . (3.3) 



For k — 1 we have that 



where /3[uj] = \/l — uj 2 . We then find that the first derivative of H p with respect to u/ evaluated at uj' = ui is indeed 
zero. The second derivative evaluated at u/ = uj leads to the following expression: 



d 2 H p 
duj 12 



(7T^(^-3)+4tanh- 1 ( v /^)) 
g 2 uj 2 {uj 2 — l) 2 



(3.5) 



This function is zero at uj = 0.697586 and the second derivative is negative below this value of uj showing an instability. 
In our numerical simulations of the unforced NLDE |17j . we find that below this value the solitary waves are metastable, 
with the time for the instability to set in depending on the value of ui. 



IV. NLDE IN EXTERNAL ELECTROMAGNETIC FIELDS 

We add electromagnetic interactions through the gauge covariant derivative 

then under the combined transformations 



1 



-> A„ - ^A, * -> y e - zA{x) 



(4.1) 



(4.2) 



the Lagrangian is invariant. Again the conserved current is given by Eq. (2.39). The gauge invariant Lagrangian for 
the external field problem is 



L 



M [^a,,* - d^-y^] - to** + _ eVyfA^ 

— J hi ~\~ _1_ 



(4.3) 



One again finds that energy-momentum is conserved, with E and P given by eqs. (2.43) and (2.44 1. Because of 



Lorentz invariance it is sufficient to calculate the energy-momentum tensor in the comoving frame v — 0. The energy 
-momentum tensor in an arbitrary frame is then given by Eq. (2.45). Using the freedom of gauge invariance, one can 
choose the axial gauge A 1 = 0, eA — V(x). In the axial gauge the Dirac equation becomes 



r/^* - to* + 5 2 (**) K * - ~f°V{x)^ = 0. 



(4.4) 



Going into the rest frame and choosing for tpo the representation of Eq. ( 2.5 ), we find that the Dirac equation becomes 



d x A + (to + ui)B - g 2 [A 2 - B 2 \ K B - V(x)B = 0, 
d x B + (m - u))A - g 2 [A 2 - B 2 \ K A + V(x)A = 0. 



(4.5) 



As in the case of the standard NLDE equation, we can again look for static solutions of the form i/jq, where 
A = Rcos8 and B — RsinO. The conservation of T 11 for static solutions that vanish at infinity yields the equations 



tuip't/j — mipip + Lj = 0, 



(4.6) 



where 



K + l 



V{x)^ip. 



Multiplying the Dirac equation [Eq. (4.4)] on the left by * and using Eq. (2.3) we obtain 



uj^ip + iip^d x ip - m^ip -I- g 2 (ipip) K+1 - V(x)ip^tp = 0. 
The energy density is given by 



(4.7) 



(4.8) 



[*7 X X * - dx^j 1 ^] + to** 



k+1 



hi + h 2 



(4.9) 
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From Eq. (4.6 1, (4.7) and Eq. (4.8) we again find that 



h 3 = -hi 



(4.10) 



Multiplying the Dirac equation [Eq. (4.4)] on the left by ^ and using Eq.(2.3), eliminating the self-interaction term 
using Eq. (??), we obtain 



iip"/ 1 d x ip + k \mij)ip + V(x)ip*ip — ujtp'tp] = 0. 



Next using the ansatz Eq. (2.5): 



'</'o 



R(x) cosd 
i R(x) s'md 



we obtain the first order differential equation for 8, namely 



(10 

dx 



-ku! + Kmcos(28) + kV(x), 



and for R(x) we obtain from Eq. (4.6) 



(4.11) 



(4.12) 



(4.13) 



[u - V(x)]R 2 - mR 2 cos 28 + -^—(R 2 cos 26) K+1 = 0, 

K + 1 



or 



R z 



(k + l)(mcos26> - u> + V(x)) 
5 2 (cos26») K + 1 



1/k 



(4.14) 



(4.15) 



Notice, if we now choose V{x) = \x cos 28 we arrive at the solutions found earlier with m — > m + /i. Note that for 
a bound state without the external potential we needed U) < m. So now we need u) < m + ji for (3^ to be real. The 
potential is then 



V(x) = fx 



(to + n) + u cosh 2/3kX 
oj + (m + (j,) cosh 2/3kX ' 



where 



@k = n\J (m + fi) 2 - lo 2 



(4.16) 



(4.17) 



The eigenvalue uj is set by fixing the charge Q of the solitary wave. We can choose a negative \i satisfying ui < m + fi 
so that the potential at small x looks much like a harmonic trap. Namely, choosing n = 1, /i = —1/4, to = 1, u> = 1/2, 
then 



2 cosh 



V(x) = - 



i 



12 cosh 



(4.18) 



This is plotted in Fig. [T] The charge density for the solitary wave corresponding to this external potential is plotted 
Fig- i 



At small x, V(x) has the expansion: 



V(x) 



1 x 2 13x 4 



O x b 



4 32 1536 
Another solution is found if we let 

V{x) = fism(26(x)). 

Then the solution to the differential equation: 

8'(x) — — Kid + Kmcos(26(x)) + K^,sin.(28(x)) 



(4.19) 

(4.20) 
(4.21) 
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FIG. 1: (Color online) Potential V(x) corresponding to Eq. ( |4.16[ ) with it = — 1/4, uj = 1/2. 




FIG. 2: (Color online) Charge density R (x) corresponding to the potential V(x) given by Eq. (4.161 with = —1/4, uj = 1/2 



is found to be 

6{x) = tan" 
Letting xq = 0, we have that 



/ /i + \J n 2 + m? — uj 2 tanh (x — xq) \J fi 2 + m 2 — uj 2 ^j 



m + uj 



tan 9 



fi + \J /j, 2 + m 2 — cj 2 tanh ^ki\/ /i 2 + m 2 — uj 2 ^ 



m + uj 



The potential is thus given by 



V 2 (x) = /isin20 = 2/z 



tan 



1 + tan 2 

2/i ^/i + -\//i 2 + m 2 — cj 2 tanh ^Kz-y/ /i 2 + m 2 — 

i ^M+a/ U 2 +m 2 — u 2 tanhffcx-v/ u 2 +m 2 — uj 2 \ \ 



(4.22) 



(4.23) 



(4.24) 



For the values K = 1, /U = — 1/4, m = 1, uj = 1/2, we find that Vz(x) has a kink-like form shown in Fig. [3] The charge 
density R 2 corresponding to these values of the parameters in V^x) is shown in Fig. [4] 
For small x for these parameters 



Tr/S 3 1365x 54249a; 2 „ , , N 

V^(ac) ~ h O (x 3 ) , 

v y 37 5476 810448 v 7 ' 



(4.25) 



if instead we choose /i = 1/4 and leave the other parameters the same we get the opposite type of kink shown in Fig. 
Now for small x we have 
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FIG. 3: (Color online) Potential V^ix) corresponding to the exact solution Eq. (4.241 with uj — 1/2, /j, = —1/4 




FIG. 4: (Color online) The charge density for the solitary wave in the potential V-i(x) given Eq. (4.24| with to = 1/2, /i = —1/4 



Tr . . 3 1365a; 54249a; 2 „ , , N 

V 2 (x) ~ 1 h O (an . 

' 37 5476 810448 v ; 

We can also solve for the more general case 

V(x) = m cos(26>(a;)) + ^ 2 sin(20(x)), 
then the solution to the differential equation is 

9'(x) = —klj + «(m + fi%) cos(28(x)) + kh% sin(20(x)), 



(4.26) 

(4.27) 
(4.28) 



vw 





0.2 






0.1 




-2 


-1 / 


1 2 




/ -0.1 





FIG. 5: (Color online) Potential V(x) corresponding to the exact solution Eq. (4.24 1 with fi = 1/4 
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FIG. 6: (Color online) Potential Vs(x) corresponding to Eq. (4.301 with uj = 1/2, /ii = /ta = —1/4 



so that we just need to change to — > to + fii, and /i —¥ H2 to obtain 

/i 2 + \J h\ + (m + Hi) 2 - uj 2 tanh ^/cx-^//^ + (to + Hi) 2 — uj 2 ^j 



tan 8 = 

The potential is then given by 



m + Hi + uj 



Vs(x) = /L«i cos 20 + /z 2 sin 28 = Hi 
N{x) 



1 - tan 2 ( 



2H2 



tan9 



1 + tan 2 8 1 + tan 2 



where 



iV(x) 



2^2 + V A*i + (Mi + m ) ; 



w 2 tanh 



(kx^J h\ + {hi + m) 2 - uj 2 ^ J 



/ii + m + cj 

/ f /i2 + Va*2 + (a*i + rn ) 2 ~~ ^ 2 tanh ^kx-^//^ + (Mi + m ) 2 ~~ <^ 2 ) J 

(/Zi + m + uj) 2 



and 



£>(*) 



[h2 + Vm 2 + G"i + m ) : 



w 2 tanh 



/Lt| + (/i! + to) '■ 



1. 



(4.29) 



(4.30) 



(4.31) 



(4.32) 



(/ii + m + uj) 2 

For w = 1/2, Hi = A*2 = —1/4, ft = 1, m = 1, V^x) has the form of Fig. [6j The charge density corresponding to this 
potential is given in Fig. [7J 

One can also consider the potential V4 — /itan 2 8. In that case on substituting t = tan0 we obtain the differential 
equation for t 



dt 

dx 



k (t 2 (fi — to — uj) + to + nt — 



(4.33) 



Thus we obtain 



KX 



s/2-s/h tan 1 



\J — (/i+Lj) 2 +m 2 +2m(a;— 3^) — m— u 



y/ (h + UJ ) 2 + 172,2 + 2m(uj — 3(i)\J H ~ \Am + w ) 2 + ™ 2 + 2m(w — 3/i) — m — uj 



\p2J~H tan 1 



^/(/i + w) 2 + m 2 + 2m(w — 3/i) v/// + \J (h + w) 2 + to 2 + 2m(w — 3/i) — to — uj 



(4.34) 
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FIG. 7: (Color online) Charge density for the solitary wave in the potential V3 (a;) corresponding to Eq. (4.301 with lo = 
l/2,/xi=M2 = -l/4. 



A. Stability of exact solutions with k = 1 in the presence of an external field 

Again we follow the method of Bogolubsky [TE\ and consider stability to changes in frequency keeping the charge 
Q fixed. That is if we parametrize a rest frame solitary wave solution by 



^ s {x,t)= Xs {x,u))e-^\ 
then we choose our perturbed wave function to be 



(4.35) 



i>[x,t,u}',tJ\ 



, X s{x,J)e- iut 



(4.36) 



Again inserting this wave function into the Hamiltonian we get a new probe Hamiltonian H p depending on both u/ ,u>. 
We will consider the case where the potential has the form 



V(x) — Hi cos 20 + fi2 sin 20. 



The energy density is given by the four terms in Eq. (4.9|, (hi + h 2 — h% + /14), where 



dO 

hi =R 2 — ; h 2 = mR 2 cos26, 
ax 

h 3 = |( J R 2 cos26») 2 ; h 4 = R 2 (m cos 26 + n 2 sin 20). 

We also have that h^ = hx. Following our previous discussion, the probe Hamiltonian has the form: 

H p [w'M = Hi[J] (/(o/.w) 2 - jV,w) 4 ) + (H 2 + H 4 )[u']f(uj',u J ) 2 . 
For our choice of potential, we have the following relationship: 

(16 



dx 



Solving this equation we obtain 



(to + fii) cos 20 + fi 2 sin 20 — w. 



T = tan = D + E tanh Bx, 



where 



We also have 



B 



^J{m + hi) 2 + ii 2 -u 2 - E 



B n ^ 2 a , , 

-r\ D = — ; di=m + fii+uj. 

di di 



2T 1-T 2 
sin20=- — T ; cos 20 = 



(4.37) 



1 + T 2 ' 



1 + T 2 ' 



(4.38) 
(4.39) 
(4.40) 
(4.41) 

(4.42) 
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therefore 



The charge is given by 



thus 



Now since 



and 



we find 



and 



where 



r 

2B 2 



(m + fix) cos 29 + H2 sin 26 — lu 
cos 2 20 



Q[u] 



2 Jl + iD + Eta^BxTl 

dig2 secn ax ^ _ (D + E tanh Bx yY ■ 



2 _ 2^ (m + pti) 2 + /j 2 , -"cj 2 " 

— 2f,.,2 _ ,,2^1 



5 2 (w 2 - Ml) 



/ii = i? 2 [(m + /ii) cos 26* + sin 29 — u] 



h 2 + hi — R 2 [(m + /ii ) cos 20 + ^ 2 sin 20] 



Hi = -uQ + H 2 + H 4 



H 2 + Hi = (m + Hx)Ti + fi 2 T 2 , 



Ti = J dx R 2 cos 20 = 2(tanh~ 1 (£ - D) + tanh _1 (.E + D)) 

2 _i / + + m2 + 2 ^im - w 2 \ 

= tanh ■ , 

9 \ Mi + m / 



(4.43) 

(4.44) 

(4.45) 

(4.46) 
(4.47) 
(4.48) 
(4.49) 



(4.50) 



Thus 



To = dx R 2 sin 20 = 



2^ 2 (m + ^!) 2/x; 



tanh 



g 2 ^ 2 -^ 2 ) ff 2 V(m + Ml ) 2 + (Mi-^ 2 ) 



flp[w'] = -a/Q[u/]/(u/, c) 2 (l - /(a/, c) 2 ) + /(«', oj) 2 (H 2 [lu'} + Hi[co'})(2 - jV, w ) 2 ). 



(4.51) 



(4.52) 



If we take the first derivative of this probe Hamiltonian with respect to u/ and set w' = w we find that this is zero 
only when fi 2 = 0. Thus the probe Hamiltonian is stationary to lu variations only when fi 2 = 0. However, we can use 
this method for the external potential V\{x). In that case the first derivative is automatically zero when lu' = lu. The 
second derivative changes sign at a particular value of lu — lu* , which for to = 1 is the solution to the equation: 



2(Ml 



lrcoth" 



Mi + l 



vV? + 2^ - lu 2 + 1 



3M? 



6mi 



3 = 0. 



(4.53) 



Now we have that for there to be an allowed real solution lu < to + /i. So if fi is positive, the original lu c is shifted 
upward and so is the region in lu space where the solutions are real. If /i is negative lu c is shifted downward, but also 
the allowed regime of lu is decreased . The net result is that the possible stability of stability is approximately 30% of 
the allowed region for real solutions independent of the value of /i. However, preliminary simulations show that some 
of the exact solutions with the potential (4.16) are metastable. 
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V. VARIATIONAL ANSATZ FOR THE NLDE IN EXTERNAL FIELDS 



The gauge invariant Lagrangian for the external field problem is 

L = (^] [§7^^* - ^§7^*] - m** + -^—(^) K+1 - e^A^ . (5.1) 
y 2 J k -\- 1 

Using the freedom of gauge invariance, one can choose the axial gauge A 1 = 0, eA = V(x). 

Our ansatz for the trial variational wave function is to assume that because of the smallness of the perturbation 
the main modification to our exact solutions to the NLDE equation without an external field is that the parameters 
describing the position, momentum, boost and phase become time dependent. 

That is, we replace 

vt ->• q(t); T] ^> r](t); juv p(t); wt' = ju(t — vx) ->• (j>(t) — p(t)(x — q(t)), (5.2) 

where 4>{t) = ujjt — p(t)q(t). 

Thus our trial wave function in component form is given by: 



*i(M) = (cosh |v4(.x') + isinh|B(a;')) e-^ +ip{x - q) , 
V 2 (x,t) = (smh^A(x') + icosh^B(x')^ e ~^ +l ^ x - q \ 



where x' — coshr/(t) (x — q(t))- Using this trial wave function we can determine the effective Lagrangian 
variational parameters. Writing the Lagrangian density as 



where 



Here B'(x') 



where 



£ = £ 1 + £ 2 + £ 3 , 

A = \ (*7 M ^* - d^V) 
= Re (i¥d t ^ + i^axd^) 

= (pq + <$> ~ P(x - ?)) {A 2 (x') + B 2 {x')) cosh 77 + q sinh 77 cosh 77 (A(x')B'(x') - B{x')A'{x')) 
+ (x -q)r] sink 2 -q(B A' - AB') -psinhr/^ 2 + B 2 ) - cosh 2 77 (AB' - BA') . 

dB(x') 



(5.3) 
for the 

(5.4) 



d , . Integrating over x and changing integration variables to z = (x — q) cosh 77 one obtains 
Li = / dx£\ = Q I pq + <j) — p tanh 77 ) — /q (cosh f] — q sinh 77) , 



= / dz[A 2 (z) + B 2 (z)}, 



is as given by Eq. (2.41 ). The second integral is 

I = J dz {B'A - A'B) = H l 



where Hi is the rest frame kinetic energy and is given by Eq. (2.51) 
Now let us turn to L 2 - We have 



£■2 = -to** 
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K+l 



K+l . 



** = A 2 (x')-B' z {x'), 



thus 



Li 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



£odx 



-h 



cosh 77 (« +1) cosh 77 



(5.9) 



(5.10) 
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where 

h = J dz(A 2 (z)-B 2 (z)), (5.11) 

h = J dz(A 2 (z)-B 2 (z)) (K+1) . (5.12) 

Finally 

C 3 = -eA *7°* = -V(x)&$ = -V{x) coshr][A 2 (x') + B 2 {x')] = -V(x) cosh 77/3(2;'), (5.13) 

thus 



dzp(z)V 



z 



cosh 77 



+ q(t) 



= -U[r,(t),q(t)}. (5.14) 



Putting these terms together we obtain: 

L = Q(pq + 4> — p tanh 77) — Iq (cosh rj — q sinh 77) 



™Hi + ? THY >, ^ - U[r,(t), q(t)}. (5.15) 

cosh 77 (« + l)cosh?7 

We now get the following Lagrange's equations: 

d 5L dQ , 

--^=0^-=0^Q = const, (5.16) 

i.e. the charge is canonically conjugated to the phase <j). The canonical soliton momentum, which is conjugated to 
the soliton position, is 

SL 

P q = -77- = QP + I a sinh 77; 
5q 

dP„ T , . 5L dU 

— f = Qp + I cosh rjr] = — = - — . 
dt oq oq 



From 



SL 

— = -> q = tannr?; (5.17) 
dp 



which implies sinh 77 = jq and cosh 77 = 7 = (1 — q 2 ) 1 . 

= — > — Qp scch 77 — / (sinh 77 — <j cosh 77) 



£77 



o 2 \ 9C7 

+ tanh?7seclrn( mh - — ^ — -I 2 - tt" = °- ( 5 - 18 ) 

/ or/ 



Changing variables to 4 = tanh 77 and using 

c£t7 



^ scch 2 ??, (5.19) 



we obtain 



and 



a 2 \ dU 

■^-wvTjy-i* <5 - 20) 
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From Eq. (5.20) we also have 



Qp 



dt 



mil 



d dU 



(k + 1) ) dt dq 



-I 



Combining Eqs. (5.22 1, (5.21) we obtain: 



d(iq) „ r 

v^jf - = Feff[q,q], 



where 



fi = mil - + I 2 + h ; Fe// [q, q] 



d^dU _ dU 
dt dq dq 



(5.22) 



(5.23) 



(5.24) 



Now for the NLDE without the presence of external forces, the solitary wave in the frame with v — obeys the 
relationship 



uipty - rmpip + (iptp) K+1 = 0. 



For our problem this converts into 

u(A 2 + B 2 ) - m(A 2 - B 2 ) + t^—(A 2 - B 2 )^ = 0. 
v ' y ' (k+1) 



Integrating this relationship we obtain: 



mil 



thus we can write Eq. (5.23) as 



where 



(k + 1) 
d(Mq) 



dt 



= F eff[q,q], 



Here F e ft is given by Eq. (5.24), where 

U{q,q) 



M = (Qui + Ioh = M07. 



dz V[ q+^ ) [A 2 (z) + B 2 (z) . 



and 7 = 1/yl — q 2 . We can rewrite the relativistic force equation as 

j 3 M q = F eff [q,q]. 
Using the rest frame identities of the Appendix we have that 

M = I Q + ujQ. 



(5.25) 

(5.26) 

(5.27) 

(5.28) 
(5.29) 

(5.30) 

(5.31) 
(5.32) 



Here Jo is given by Eq. (5.8), and Q is given by Eq. (2.41). 

It is useful to rewrite the equation for the canonical momentum P q using the definition of M and Eq. (5.20) as 
follows: 



P q = Qp + Ioiq = M 7<? 



dU[q,q] 



(5.33) 
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A. Energy-momentum tensor 



The fact that the external potential is explicitly independent of time means that the energy of the solitary wave is 
independent of time. The energy density is given by 



Straightforward integration leads to 

E = [ dxT 00 = Qpq + 7 / + ~h - , g —r, h + U[q, q]. 



Using the identities Eq. (5.27) and Eq.(5.29) we can rewrite this as 

E = A/07 + Qpq - l^Qq 2 + U[q, q]. 



From (5.20) and (5.27) we have that 



Qp = 'fquQ - 



dU 
dq ' 



thus we can write the energy of the solitary wave in the convenient form: 

.dU 



E = A/07 + U[q, q]-q 



dq ' 



(5.34) 



(5.35) 



(5.36) 



(5.37) 



(5.38) 



The conservation of energy will be important to test our numerical integration schemes in Sec. VII. 

For time independent external forces the total momentum of the solitary wave is not conserved but changes de- 
pending on the external force. We have that 



Explicitly we obtain 



P = I T Q1 dx = ~ I dx(^d x ip - d x if>^). 



P = iqlo+pQ, 



where Iq is given by Eq. (5.8). Using Eq. (5.20), we can rewrite this as 

. dU[q,q] 



p = iM 



dq 



(5.39) 



(5.40) 



(5.41) 



which we recognize as identical to the canonical momentum P q = given by Eq. (5.33). The Lagrange equation for 



P q is 



p = dU[g,g] = £ / _ dU[q,q} 

dq dt \ dq 



(5.42) 



VI. STABILITY CONJECTURE 



For the NLSE it was shown earlier [91 [TT] that a reliable dynamical stability criterion for the breakup of the solitary 
wave under external forces was that the solitary wave will be stable if 



dp(t) 



> 0. 



(6.1) 



Here p(t) is the normalized momentum of the solitary wave P(t) /M(t), where M = j dx^*(x, t)^(x, t) is the "mass" 
of the solitary wave. For the NLDE Q takes the place of M . However, Q is a conserved variable so one can use the 
canonical momentum P(t) instead of P(q,q)/Q to study stability. Using our collective coordinate theory, this leads 
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to the criterion that a necessary (but not sufficient) condition for stability of the solitary waves of the CC theory is 
that 



Note that the r.h.s. of Eq. (6.2) plays the role of a time dependent mass. The sufficient condition for the solitary 

ations is that 

< 0. (6.3) 



wave solution to the CC equation to be unstable in our simulations is that 

dP(q,q) 



dq 

Following Comech's reasoning [14J . we expect that in the nonrelativistic regime where u is close to m that this 
criteria will be valid in determining stability. 

VII. SIMPLE EXTERNAL POTENTIALS 
A. Simulations 

The numerical simulations have been performed by means of a 4th order Runge-Kutta method. We choose N + 1 
points starting at n = and vanishing boundary conditions \?(±L,i) = 0. The other parameters related with the 
discretization of the system are x E [—100, 100], Ax = 0.02, At = 0.0001. For our initial conditions on the soliton, we 
use the exact 1-soliton solutions of the unforced nonlinear Dirac equation discussed in Sec. II. Since we would like to 
compare the exact numerical solution with the solution of the CC equations, we need to define how we determine the 
position of the soliton. In our numerical computation of q(t) we have used the first moment of the charge, i.e. 

q(t) = Q 1 (t)/Q(t), (7.1) 

where 

«,= /*,#.(,)•(,), «-/*,..<,).(,), (7.2) 
B. linear potential (ramp potential) 



Consider the constant external force with scalar potential V(x) = —V\x 1 and V\ > 0. We then have from Eq. (5.14) 
that 

U = -V iq {t)Q, (7.3) 

the force law then becomes 

j t (M ol q)=V 1 Q. (7.4) 

Integrating once (starting at an initial velocity q ) one has 

g 

where c\ = ViQ/Mq and c 2 = , q ° , ■ Integrating we obtain 



c 1 t + c 2 , (7.5) 



y/( Cl t + c 2 ) 2 + 1 Vc 2 2 + 1 
?(*) = - 1-9(0)- (7.6) 



Cl Cl 



This is the standard result for a relativistic point particle undergoing constant acceleration. If we choose q(0) 
q(0) = 0, we get the simpler expression 
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The energy of the solitary wave is just 



E = M07 + U = M , 



and the force law is now 



so that 



aq 



P = ViQt. 



Since U in this case is independent of q, we find from Eq. (6.2 1 that 



DP 
dq 



rM > 0, 



(7.8) 
(7.9) 
(7.10) 
(7.11) 



thus the necessary condition for stability is fulfilled. In this section and what follows we will confine ourselves to the 
case where k = 1 and also Q = 1. For that case from the Bogolubsky stability requirement |16j we know that without 
forcing when uj < uj c = 0.697586, the solitary waves are unstable. 

Let us look at the case where the unforced solitary wave is stable. For u) = 0.9, g = 1 we have solved numerically 
the NLDE for V\ = 0.01 (as well as V± = 10~ 3 and V\ = 10~ 4 ). We find for all these values of V\ the solitary wave is 
stable at all simulation times and t he center of the solitary wave follows the analytic formula we derived from the CC 
equation for q(t), namely Eq. ( 7.7 1. In Fig. [8]we display the results of the simulation for p(t),q(t),P(t) for V\ = 0.01. 
We notice that the width of the soliton gets Lorentz contracted as the velocity increases (this effect is not apparent 
for the smaller values of V\). Because the charge is conserved, the height of the solitary wave increases due to the 
increase of j(t). 

For the case uj = 0.3, the unforced solitary wave has double humped behavior and is unstable at late times. Here 
our simulations show that until the instability sets in (around t s» 110, for m = l,g = 1,V\ = 0.01) the position 
of the solitary wave follows the analytic solution of the CC equation Eq. (7.7 1. However, the actual shape of the 



solitary wave becomes asymmetric with the left hump becoming higher than the right hump as a precursor for the 
wave becoming unstable. This is shown in Fig. |9j where PQ(x,tfi X ed) is plotted against x for various t — t*. In Fig. 
10 we give results of the simulation for the case where uj = 0.3, V\ = 0.0001. Here, looking at q(t) we explicitly see 
that around t — 120, the solution of the NLDE diverges from the solution of the CC equation. Also for this value 
of the potential the solitary wave humps are symmetric and that the single solitary wave breaks up into two solitary 
waves with some radiation when it goes unstable. 



C. Harmonic Potential 



Let us consider the case of an external harmonic potential, V{x) = \V^x , V2 > 0. For that case from Eq. (5.14) 



we find that 



U =-Vo dz 



q+-) [A 2 {z) + B\z) 
7 



-V 2 [Qq 2 {t) + {l-q 2 {t))h] , 



where I 3 = J^dz z 2 [A 2 (z) + B 2 (z)}. From Eq. (5.24) we have 



Feff[q,q\ = -V 2 hq-V 2 Qq 



leading to the equation of motion [see Eq.( 5.28) 

d 



dt 



{[M o7 [«(*)] + V2I3] q} = -V 2 q(t)Q. 



(7.12) 



(7.13) 



(7.14) 



This can be rewritten as 



[AW + V 2 I 3 ]q + V 2 Qq = 0. 



(7.15) 
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FIG. 8: Simulations with a ramp potential, V(x) — —V\x. Upper panel: Charge density at t* = 16.6; 33.3; 50; 66.6; 83.3; 100. 
Middle and lower panels: soliton position q(t) and momentum P(t) from analytical results of the CC equations (solid lines) and 
from numerical simulations (dashed lines) of the forced NLDE. The curves are super-imposed. For the final time of integration 
the relative error of q(t) is of order 10~ 5 . Parameters: g = 1, m = 1, u> = 0.9 and Vi = 0.01. Initial condition: exact soliton of 
the unperturbed NLDE with zero initial velocity. 
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FIG. 9: Simulations with a ramp potential, V(x) — —V\x with ui in the unstable regime. Left upper panel: Charge density at 
t* = 16.6; 33.3; 50; 66.6; 83.3; 100. Left lower panel: Charge density at t* = 110. Right panels: q{t) and Pit) from analytical 
results of the CC equations (solid lines) and from numerical simulations of the forced NLDE (dashed lines), the curves are 
super-imposed. Parameters: g = 1, m = 1, w = 0.3 and Vi = 0.01. Initial condition: exact soliton of the unperturbed NLDE 
with zero initial velocity. 



In the nonrelativistic regime where 7[(j(£)] w 1 we recover the oscillator equation for the collective coordinate q(t) 
namely 



q + n 2 q = o, n 2 



M + Vih 

Note that the rest mass is increased by the term V2/3 > 0. For initial conditions qo = 0, q(0) = v Q we obtain 

q(t) = ~ sinfii. 



(7.16) 



(7.17) 



1. Energy conservation 



From the energy conservation equation in the form of Eq. (5.38) 



r , dU 

E = M ol + U[q,q]-q — 



(7.18) 
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FIG. 10: Simulations with a ramp potential, V(x) = —V\x with uj in the unstable regime. Upper panel: Charge density at 
t* — 0; 50; 100. Middle panel: Charge density at t* = 150. Lower panel: q(t) from analytical results of the CC equations (solid 
line) and from numerical simulations (dashed line) of the forced NLDE. Parameters: g = 1, m = 1, u = 0.3 and Vi = 0.0001. 
Initial condition: exact soliton of the unperturbed NLDE with zero initial velocity. 



we obtain 

E = Moj + \QV 2 q 2 + \v 2 h{l + q 2 ). (7.19) 

In the low velocity limit, we need to keep the first two terms in the expansion of 7 in the expression for the energy; 
namely (here we suppress the speed of light in the v /c expansion) 

7 = 1 + i<j 2 + .... (7.20) 
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We then have for the solution q(t) — j£ sin fit the usual equipartition of energy and that the non-relativistic energy 
is twice the initial kinetic energy (apart from a constant) 

E = M + 2T; T = {M + V 2 I 3 )v 2 (7.21) 

with the effective mass of the solitary wave increased over the unforced case by the quantity V2/3. 



2. Canonical momentum and stability Criterion 

From the equation for the canonical momentum, Eq. ( 5.41[ ), we find 

P = (M + V 2 I 3 )q = (M o7 + V 2 h)q. 



This again shows the mass increased by V2/3. The stability criterion Eq. (6.2 1, leads to 

9P{qA) a M , v j . n 

— = <-y Mo + V2-L3 > U. 



(7.22) 



(7.23) 



Thus the necessary condition for stability of the solitary wave is fulfilled. 

We would now like to see how well the CC equations for q and P, namely Eqs. ( |7.15 ) and (7.22 1, compare with the 
numerical solutions of the forced nonlinear Dirac equation. We will choose our initial condition to be q(Q) = 0, q(0) = v 
and study both the nonrelativistic regime (vq — 0.1) and the relativistic regime (vq = 0.9). The external potential 
can be written as 



V = 2 V " = 2 



(7.24) 



which identifies the characteristic length of the potential as I — l/yT^. We would like to choose the characteristic 
length of the potential to be large compared to the width of the solitary wave which is 1/(2/3), (3 = y/l -cj 2 . Choosing 
V-2 = 10 -4 accomplishes this requirement. At low velocities both P(t) and q(t) are proportional to cos fit, thus P{q) 
is a straight line with positive slope. 

First let us consider the regime where the unforced solitary waves are stable, and choose lu = 0.9, g = 1. In the non- 
relativistic regime (vq — 0.1) we get the results shown in Fig. 11 The oscillations of both q(t) and P(t) are harmonic 
as predicted by the CC equations. The charge density maintains its shape as its position oscillates periodically in 
time. For vo = 0.9, q and P again follow the CC equations for a little less than half the oscillation period but then 
the solitary wave becomes unstable and the exact simulation of q and P then diverges from the solution to the CC 
equations. This is shown in Fig. [12] 

Next we consider the regime < lu < lo c — 0.697586 where the unforced soliton is unstable. For the parameters 
to = 0.3, g = I.V2 = 10~ 4 , we obtain the typical result found in the unforced problem that at around t = 120, the 
solitary wave becomes unstable. Until then the CC equations for q{t) and Pit) track well the exact solution. This is 
seen in Fig. |13| However, the wave function starts becoming asymmetric at late times and departs from our symmetric 
ansatz even before the soliton becomes unstable and breaks into two solitary waves plus some radiation. 
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FIG. 11: Harmonic potential, V(x) = (V 2 /2)x 2 . Upper panel: Charge density at t* = 0; 133.3; 266.6; 400; 533.3; 666.6; 800. 
Middle and lower panels: soliton position q{t) and momentum P(t), from the numerical solutions of the CC equations (solid 
lines) and from numerical simulations (dashed lines) of the forced NLDE. The curves are super-imposed. The charge Q — 0.96864 
and energy E = 0.93921 are both conserved. Parameters: g = 1, m = 1, u = 0.9 and Vi = 10 -4 . Initial condition: exact 
soliton of the unperturbed NLDE with initial velocity v(0) — 0.1. 
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FIG. 12: Harmonic potential, V(x) = (V2/2)x . Relativistic regime, unstable soliton. Upper panel: Charge density at 
t* — 0; 133.3; 500. Middle and lower panels: soliton position q(t) and momentum P(t), from the numerical solutions of the CC 
equations (solid lines) and from numerical simulations (dashed lines) of the forced NLDE. The curves in P(t) are super-imposed 
only till t = 370. Parameters: g = 1, m = 1, w = 0.9 and V2 = 10~ 4 . Initial condition: exact soliton of the unperturbed NLDE 
with initial velocity v(0) = 0.9. 
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FIG. 13: (Color online). Harmonic potential, V(x) — (V2/2)x 2 with u in the unstable regime. Left upper panel: Charge 
density at t* = 0; 133.3. Right upper panel: Charge density at t* = 150. Lower panels: soliton position q(t) and momentum 
P(t), from the numerical solutions of the CC equations (black solid lines) and from numerical simulations (dashed lines) of the 
forced NLDE. The energy (red or middle curve) and charge (blue or upper curve) are also plotted. Parameters: g = 1, m = 1, 
lj — 0.3 and Vi = 10 -4 . Initial condition: exact soliton of the unperturbed NLDE with initial velocity v(0) = 0.1. 



D. Spatially periodic potentials 

Next consider a spatially periodic potential 

V(x) — — ecos kx, e > 0, 



(7.25) 



where the spatial period L = 2ir/k ^> 1//3, and 1/(3 is the width b of the solitary wave. The potential is then a 
function of q, q and is given by 



/kz 
dz cos — [^4 2 [z] + -B 2 [z]] = — e cos kq I&[q\ 



so that 



dU 
dq 



ke sin kql^ [q] , 



(7.26) 



(7.27) 



dU 
dq 



/kz 
dz zs'm — [^4 2 [z] +S 2 [z]] = — kejq cos kq I§[q\. 



(7.28) 



2N 



The generalized force is obtained from 



„ r ., d dU dU 



We find that 



where 



d^dU 
dt dq 



T = -ke cos kq q [j 3 h[q] ~ k(-/q) 2 1 6 [q\] 



h[q] = J dzz 2 cos ^ [A 2 [z] + B 2 [z}] . 
The generalized force equation can be written as 

F[q, q] = -keI A (q) sin(kq(t)) - ke cos kq q [j 3 I 5 [q] - Hiqfh [q\] 
In the nonrelativistic limit (q) 2 < 1, 7 w 1 and we obtain for the force law: 



= M -( 7 g)- 



where 



(Mq + feelg cos kq)q + kel® sin kq = 0, 



ll = j dz cos kz [A 2 [z] + B 2 [z\] ; 1% = J dz zsinkz [A 2 [z] + B 2 [z\] . 



(7.29) 



(7.30) 



(7.31) 



(7.32) 



(7.33) 



(7.34) 



When the potential is weak (e <C 1) then kel® <C Mo and we obtain the pendulum (or sine-Gordon) equation 

M q + kell sin kq = 0. (7.35) 
Letting C = kel® /Mq , the solutions are given by 



2 /l 



q(t) = T am -VkJ 2Ct 2 + kc^ 2 + ACc 2 t + 2kc Y c 2 t + 2Cc 2 + fccic 



AC 



(7.36) 



k V 2 "" v ~" ' '"" x " ' ' ' ~ 2 ' 2C + fcci 

where ci, C2 are integration constants to be determined by the initial conditions, g(0) = 0; q(0) = v . Here 

am[«, i] = JacobiAmplitude[M, I], (7.37) 
where the modulus parameter I (usually denoted by m) is I = 2 c 4 ffc Cl ■ For the above initial conditions we find: 



q(t) = -am 



kv Q t AC 
~ 1 T , ~k^ 2 



q(t) = — sm sn 
k 



kv t AC 



(7.38) 



J. Energy conservation 



From Eq. (5.38) we have that the soliton energy is given by 

E = 7M0 — e cos [<j] + e cos fcg kjq 2 I$ [q] . 



In the nonrelativistic limit we obtain 



, / M o 



E = M - e cos fc<?/' 



In the case of a weak potential (except of Mq — ► for lu — > 1) 



efc/5 cos /cq^ g 2 



E = ( 1 + y ) M - e cos kql%. 



(7.39) 



(7.40) 



(7.41) 
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2. Soliton momentum and dynamical stability 



The soliton momentum is given by Eq. (5.41 ) and becomes 



In the nonrelativistic regime we obtain 



The necessary condition for stability 



P = 7 (M + kel 5 [q] cos kq) q. 



P = (M + hell cos kq) q. 



dP 
dq 



- = (M + hell cos kq) > 



(7.42) 
(7.43) 

(7.44) 



is satisfied except in the regime where Mq — > which is when uj — > 1. In that regime the soliton is very broad and 
the condition 2ir/k 3> b is not fulfilled. 



3. Numerical results for q(t) and P[t) 

For the pendulum equation there is a critical initial velocity at which the coordinate q(t) makes a transition from 
periodic motion to unbounded motion. This occurs when the modulus parameter 1 = 1. This yields the condition 



'4J°e 
M " 



(7.45) 



Depending on our choice of parameters, for small enough e, v c will be in the nonrelativistic regime. We choose oj 
to be in the stability region for the unforced problem (see Sec. III). For g = 1, m = 1, oj c = 0.697586. Choosing 
uj = 0.9, then the width of the solitary w ave is 1/(2 (3) — 1.15. If we choose k = 0.1, then the characteristic wave 
length 2n/k = 62.8 > 1/(2/3). From Eqs. ( |2.55[ ) and ( |2.56[ ) we have that 



Q = 0.968644; H x = 0.0625108 = I , 



(7.46) 



M = H x +ojQ = 0.934291. 



(7.47) 



The other constants for this initial condition from Eq. (7.34 1 are 



II = 0.94632; 7~£ = 0.433477. 



(7.48) 



We have first compared the analytical solution Eq. (7.381 of the pendulum equation with the numerical solution of 
Eq. (7.33). For e < 0.1 the results are practically identical, for e > 1 deviations occur. 

Specifically we have chosen the initial condition q(0) = 0, q = vq for the three cases (1) vq -C v c <C 1 and then vq 
slightly below (2) and above (3) the critical value v c , namely 



Vq = V c =F -001. 



(7.49) 



Choosing e = 0.001 yields v c = 0.0636619 which is in the nonrelativistic regime, so we expect Eq. (7.38) to hold. In 



Fig. 14 we show that for vq = 0.01 the analytic nonrelativistic result and the numerical solution of Eq. (7.33) give the 



same results as the solution of the NLDE. We also see that the shape of the charge density does not change in time. 
In Figs. 15 and 16 we show that just below and above the critical velocity, respectively, the analytical result (7.38) 



agrees with the numerical solution of Eq. (7.33), but both results differ very slightly from the simulation results. 



VIII. CONCLUSIONS 



We presented analytical and numerical results for the forced NLDE for the scalar-scalar interaction and arbitrary 
nonlinearity; for numerics we used k = 1. For the forcing terms we used simple test potentials such as ramp, harmonic 
and periodic potentials. We discussed one criterion for the stability of the solitons. We have given two sufficiency 



30 





0.01 



0.005 



-0.005 



-0.01, 




2000 3000 
t 



4000 



FIG. 14: Periodic potential, V(x) = — ecos(kx), very low initial velocity. Upper panel: Charge density at t* = 0; 2666.6. 
Middle panel: soliton position q(t), from the numerical solution of Eq. (7.331 (solid line), approximate analytical expression 
(7.381 (dotted line), and from numerical simulations (dashed line) of the forced NLDE. The three curves are super-imposed. 



Lower panel: momentum P(t), from the numerical solutions of Eq. (7.331 (solid line) and from numerical simulations (dashed 
line) of the forced NLDE. The curves are super-imposed. Parameters: g = 1, m = 1, u = 0.9, e 
condition: exact soliton of the unperturbed NLDE with initial velocity v(0) = 0.01. 



0.001 and k = 0.1. Initial 



conditions for instability which complement the results given in [14j which are valid only in the nonrelativistic regime 
where u> < m. For one solution the criterion is that the soliton minimizes its energy for fixed charge as a function of 
lo. Another criterion is a negative slope in dp/dq. Stability was studied numerically for various cases. 
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FIG. 15: Spatially periodic potential, V(x) = — ecos(fca;), initial velocity just below v c . Upper panel: Charge density at 
t* — 0; 5000. Middle panel: soliton position q(t), from the numerical solutions of Eq. ( 7.33 1 (solid line), approximate analytical 
expression (7.381 (dotted line), and from numerical simulations (dashed line) of the forced NLDE. Solid and dotted lines are 
super- imposed. Lower panel: momentum P(t), from the numerical solution of Eq 
simulations (dashed line) of the forced NLDE. Parameters: g = 1, m = 1, w = 0.9, 
exact soliton of the unperturbed NLDE with initial velocity v(0) = 0.0626619. 



(7.331 (solid line) and from numerical 
0.1. 



0.001 and k 



Initial condition: 



We also developed a variational ansatz for the NLDE in external fields and compared approximate analytical 
solutions with numerical solutions of the collective coordinate equations and with simulation results. In all cases, we 
found very good agreement between the three. 
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FIG. 16: Spatially periodic potential, V(x) = — €Cos(kx) y initial velocity just above v c . Upper panel: Charge density at 
t* — 0; 4000. Middle panel: soliton position q(t), from the numerical solutions of Eq. ( 7.33 1 (solid line), approximate analytical 
expression (7.381 (dotted line), and from numerical simulations (dashed line) of the forced NLDE. Solid and dotted lines are 
super- imposed. Lower panel: momentum P(t), from the numerical solution of Eq 
simulations (dashed line) of the forced NLDE. Parameters: g = 1, m = 1, u = 0.9, 
exact soliton of the unperturbed NLDE with initial velocity v(0) = 0.0646619. 



(7.331 (solid line) and from numerical 
0.1. 



0.001 and k 



Initial condition: 



The simulations in this paper are confined to the k = 1 case. The numerical stability of solitons for general k will 
be presented in a subsequent publication [17] . The semiclassical reduction of NLDE to NLSE and implications for 
soliton stability has been recently discussed in a rigorous fashion by Comech |14j . Our numerical findings agree with 
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his analysis in the non-relativistic regime. We intend to study the forced NLDE with vector-vector interactions and 
arbitrary nonlinearity in a separate paper. 
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Appendix: Rest Frame Identities 

In the rest frame, energy-momentum conservation leads to identities among the various integrals that arise. The 
Lagrangian in the axial gauge is given by 

2 

L = \ [*7% - d M *7 M *] - to* + (**) K+1 - V(x)¥$> . (9.1) 
In the rest frame, the wave function is given by: 

_ ( A(x) \ n -i u t 



* = V,e-**= ( )e- wt . (9.2) 



The energy-momentum conservation is given by Eq. (2.7) and leads to two independent equations. The first is 



d T 0x + 8 X T XX = 0. (9.3) 

In the rest frame T Qx — 0, so that T xx — constant. If the solution goes to zero at infinity then the constant is zero. 
We have then the relationship: 

T xx = % - [f 7*d a: # - d x Vj x y] + L 

a 2 - 

= unity - miptp + (ijnl)) k+1 - U(xW tp = 0. (9.4) 

k + 1 



ojQ - rnh + -^—h - I dxp(x)V(x) = 0. (9.5) 



Integrating over space we get the relations: 

K + 1 

The second conservation law is 

<9 T 00 + d x T x0 = 0, (9.6) 
which leads to the conservation of energy. The energy of the solitary wave in the rest frame defines the rest mass M 

E = J T oa dx = M . (9.7) 

We have that 

T 00 = - l -[^ x d x ^-d x ^ x ^\ 

+rmpip — (ipip) k+1 + V(x)^ip 

k + 1 

= (AB X - BA X ) + m(A 2 - B 2 ) - j^(A 2 - B 2 ) K+1 + (A 2 + B 2 )V(x). (9.8) 

k + 1 
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Integrating we obtain 



M = Iq + mil 



9 



k + r 



h + / P (x)V{x). 



Using the identity of Eq. (9.5), we then have even in the presence of interactions that 

M =I + ojQ. 



(9.9) 



(9.10) 
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